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SUMMARY 


A  data  set  which  contains  nearly  500  arrays  representing  the  measurements  of  visible  spectrum  sky  and  terrain  radi¬ 
ance  values  first  discussed  in  AFGL-TR-81-0275  has  been  organized  into  a  readily  transportable,  machine  readable  library. 
A  complementary  data  set  including  contemporary  measurements  of  atmospheric  volume  scattering  coefficients  plus 
selected  meteorological  variables  as  a  function  of  altitude  has  been  included  in  a  similar  machine  oriented  format. 

The  data  base  contents  are  described  as  are  the  storage  tape  formats  and  sample  extraction  procedures. 
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A  TRANSPORTABLE,  MACHINE  ORIENTED  LIBRARY  OF 
EUROPEAN  SKY  AND  TERRAIN  RADIANCE  DISTRIBUTIONS 
WITH  CONTEMPORARY  RADIOMETRIC  4  METEOROLOGICAL  PROFILES 

Richard  W.  Johnson 
Miriam  K.  Oleinik 


1.0  INTRODUCTION 

A  substantial  body  of  meteorological  and  visible 
spectrum  radiometric  data  describing  the  optical  state  of 
the  atmosphere  has  been  collected  and  reported  upon  by 
the  Visibility  Laboratory  of  the  University  of  California, 
San  Diego  in  cooperation  with,  and  under  the  sponsorship 
of  the  Air  Force  Geophysics  Laboratory,  Hanscom  AFB, 
Massachusetts.  Many  of  these  data  are  related  to  the 
European  environment,  and  are  available  in  hard  copy  via 
one  or  more  of  the  reports  listed  in  Table  I  I,  which  has 


TaM«  1.1  OPAQUE  related  aircraft  data  reports. 


Report  No. 

Determent 

Due  tree 

Serin 

Dale 

apoi  ML-v-acm 
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rv  v 

Review  of  Op  Poop  A  AeroaoA 

been  abstracted  from  AFGL-TR-82-0049,  Johnson  and 
Fitch  (1981).  The  OPAQUE  measurement  program  K)Pt- 
ical  Atmospheric  Quantities  in  Europe)  was  conducted 
under  the  auspices  of  the  NATO  Research  Study  Group  8 
of  Panel  IV,  AC243 

Whereas  alt  but  the  last  two  reports  listed  in 
Table  1.1  concentrated  upon  the  documentation  of  the 
vertical  structure  of  the  atmospheric  volume  scattering 
coefficient,  it  is  important  to  note  that  there  is  a  substan¬ 
tial  and  valuable  complementary  data  set  that  is  related  to 
the  atmospheric  radiance  distribution.  This  complemen¬ 
tary  data  set,  which  is  discussed  procedural^  in  Johnson 
(1981)  and  analytically  in  Johnson  and  Hering  (1981), 
contains  angularly  precise  spectral  measurements  of  the 
surrounding  4rr  radiance  Held  as  seen  by  an  instrumented 
aircraft  during  its  traverse  of  several  different  flight  tracks 
between  the  surface  and  altitudes  up  to  6  kilometers.  The 
possession  of  these  radiance  data  in  a  calibrated  and 
readily  retrievable  format  is  what  enables  one  to  character¬ 
ize  the  broad  variety  of  environmental  conditions  extant 
during  the  flight  episodes  and  thus  leads  to  the  develop¬ 
ment  of  operationally  useful  predictive  models  as  in  Her¬ 
ing  (1981). 

These  radiance  data,  in  conjunction  with  their  com¬ 
panion  scattering  coefficient  data,  are  readily  applicable  to 
the  determination  of  slant  path  contrast  transmittances, 


F1(.  1-1.  Orsuiiutionil  block  theorem 
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Table  1.1.  Summary  of  sky  «nd  terrain  nditnce  measurements. 
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T«W«  1.3.  Summary  of  aimospheric  vertical  profile  measurements. 
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atmospheric  optical  depths,  aerosol  directional  scattering 
characteristics,  flux  divergences  and  their  attendant  deter¬ 
minations  of  turbid  atmosphere,  single  scattering  albedos. 


The  description  of  a  transportable,  machine  oriented 
data  library,  shown  schematically  in  Fig.  1*1,  containing 
the  sky  and  terrain  radiances  summarized  in  Table  12, 
abstracted  from  Johnson  (1981),  and  the  selected  con* 
current  profile  data  summarized  in  Table  1.3,  is  the  sub¬ 
ject  of  this  report. 

It  is  important  that  users  of  these  library  data  take 
full  advantage  of  the  supplementary  detailed  informatioi. 
and  interpretive  guides  that  are  included  in  the  listings  of 
Table  1.1. 


2.0  DESCRIPTION  OF  DATA  BASE  CONTENT 

The  tape-oriented  data  base  that  is  discussed  in  this 
report  contains  two  independent  but  related  sets  of  meas¬ 
urements  which  are  provided  on  several  separate  storage 
tapes  whose  format  will  be  discussed  in  the  following  sec¬ 
tion  In  this  section,  we  will  confine  ourselves  to  a  brief 
discussion  of  the  pertinent  characteristics  of  the  measured 
data  itself  rather  than  that  of  the  storage  medium. 

The  two  data  seta  comprising  this  data  base  are 
described  herein  as;  a)  the  Sky  and  Terrain  Radiance 
Data,  and  b)  the  Scattering  Coefficient  and  Related 
Meteorological  Data.  For  convenience,  the  Sky  and  Ter¬ 
rain  Radiance  Data  will  often  be  referred  to  as  the 
" Scanned  data,  and  the  Scattering  Coefficient  and  Related 


P 
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Meteorological  Data  as  the  "Profiler  data.  As  a  general 
rule,  the  data  is  sorted  into  sub-sets  identified  by  flight 
number  as  indicated  in  Tables  1.2  and  1.3. 

As  with  all  of  the  radiometric  data  associated  with 
the  reports  listed  in  Table  1.1,  the  radiometric  measure¬ 
ments  contained  in  this  data  library  may  be  characterized 
by  the  illustrative  information  in  Fig.  2-1  and  Table  2.1. 
More  detailed  information  regarding  the  radiometric  cali¬ 
bration  procedures  applied  to  these  data  is  available  in 
several  preceding  reports,  AFCRL-70-0137,  AFCRL-72- 
0461,  and  AFGL-TR-80-0207  among  others  and  thus  will 
not  be  included  here. 

Users  should  also  be  aware  that  all  flight  altitudes 
associated  with  these  data  sets  are  listed  in  meters  above 
ground  level,  m(AGL),  and  have  therefore  been  corrected 
for  local  meteorological  conditions  extant  at  the  time  of 
the  measurements. 

2.1  Sky  and  Terrain  Radiance  Data 

As  discussed  in  Johnson  (1981)  and  Johnson  and 
Henng  (1981),  both  of  which  should  be  considered  com¬ 
panion  reports  to  this  current  effort,  a  data  set  containing 
nearly  300  arrays  represerting  measurements  of  sky  and 
terrain  radiance  values  has  been  developed  for  general 
application  to  image  propagation  studies  within  the  lower 
troposphere.  The  arrays  have  been  organized  for  storage 
by  flight  number  primarily,  although  date  and  site  are  also 
listed  in  the  header  information  as  will  be  noted  in  Sec¬ 
tion  3.  As  a  general  rule,  the  scanner  data  arrays  ate 
organized  as  illustrated  in  Table  2.2. 

Each  of  the  individual  arrays  indicated  <n  Table  2.2 
contains  1080  data  points  representing  either  an  upper  or 
lower  hemisphere  radiance  distribution.  Each  data  point 
represents  the  apparent  radiance  of  the  sky  or  terrain  as 
"seen"  by  the  5°FOV  of  one  of  the  airborne  scanning 
radiometers  described  in  Johnson  (1981).  There  is  a  data 
point  for  every  6°  in  azimuth  and  every  5®  in  zenith  angle 
within  each  hemisphere.  All  azimuth  notations  are  taken 
with  respect  to  the  sun  as  illustrated  in  Fig.  2-2  It  should 
be  noted  that  for  each  entry  in  the  data  summary  listed  in 
Table  1 .2,  there  are  two  data  arrays  in  the  library,  one  for 
the  upper  hemisphere  radiances,  and  one  for  the  simul¬ 
taneously  measured  lower  hemisphere  radiances. 

The  radiance  values  included  in  this  data  library  are 
the  same  as  those  used  to  create  the  graphical  displays 
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Table  2.2.  Sky  &  Terrain  Radiance  Data 
Organization  Pet  Flight 

I.  Flight  Identification.  Flight  No.,  Date,  Site,  etc. 

A.  Altitude  Number  One  (See  Table  1.2) 

1.  Spectral  Band  No.  2  (see  Fig.  2-1) 

a.  Upper  Hemisphere  Radiance  Array 

b.  Lower  Hemisphere  Radiance  Array 

2.  Spectral  Band  No.  3 

a.  Upper  Hemisphere  Radiance  Array 

b.  Lower  Hemisphere  Radiance  Array 

3.  Spectral  Band  No.  4 

a  Upper  Hemisphere  Radiance  Array 
b.  Lower  Hemisphere  Radiance  Array 

4.  Spectral  Band  No.  5 

a  Upper  Hemisphere  Radiance  Array 
b.  Lower  Hemisphere  Radiance  Array 


B.  Altitude  Number  Two 


C.  Altitude  Number  Three 

I.  ' 

j-  •  same  as  in  ’  A"  above 
4’  • 

D.  Altitude  Number  Four 


same  as  in'  A'  above 


Table  2.1.  Spectral  characteristics  summary  for  AVIZ  filters. 


Spectral  Chtrtctriftici 


Inherent  Sun  Propenia  (Johnson  (1954)] 


fteyteith  Atmosphere  Properties  <1 5*0 


FI*.  2-2.  Sky  and  terrain  coordinate  system. 


Table  2.*.  Scattering  Coefficient  &  Related  Meteorological  Data 
Organization  Per  Flight 


I  Right  Identification:  Right  No  ,  Date,  Site,  etc. 

A.  Profile  No  1  (Ascent)  (See  Appendix  C,  Fig  C-l) 

I  Scatt  Coeff.  in  Spectral  Band  No.  2  (See  Fig.  2-1) 

2.  Ambient  Temperature 

3.  Dewpoint/Frosipomt  Temperature 

4.  Relative  Humidity  (computed) 

5.  Absolute  Pressure 

6.  Atmospheric  Density  (computed) 

7.  Real  Time 


B  Profile  No  2  (Descent) 

1.  Scan.  Coeff.  in  Spectral  Band  No  3 

2.  1 

to  >  same  as  in '  A"  above 

7.  J 


Profile  No.  3  <  Ascent) 

I  Scatt  Coeff  in  Spectral  Band  4 


?} 


same  as  in  "A”  above 


D. 


Profile  No.  4  (Descent) 

1  Scatt.  Coeff  in  Spectral  Band  5 


same  as  in’ A*  above 


included  in  Johnson  (1981),  and  Johnson  and  Hering 
(1981),  and  should  be  used  with  the  same  caveats. 

THE  DATA  ARRAYS  CONTAINED  IN  THIS 
LIBRARY  HAVE  NOT  BEEN  CORRECTED  FOR  THE 
NEAR  SUN  STRAY  LIGHT  EFFECTS  AS  RECOM 
MENDED  IN  SECTION  4  OF  JOHNSON  (19*1), 
AFGL-TR-81-0275. 

An  example  of  the  clear  day  sky  and  terrain  radi¬ 
ance  distributions  that  may  be  extracted  from  these  data 
arrays  is  shown  in  Appendix  A  for  flight  C-379  These 
graphical  displays  have  been  abstracted  from  Johnson  and 
Hering  (1981)  as  a  convenient  reference  for  the  reader 
The  software  required  to  create  these  displays  is  not 
included  as  part  of  this  report 

2.2  Vertical  Profile  Data 

During  each  of  the  data  flights  listed  in  Table  1.2, 
there  were  radiometric  and  meteorological  measurements 
conducted  at  fixed  altitudes  as  well  as  during  ascent  and 
descent  modes  The  data  which  were  associated  with  the 
ascents  and  descents  of  the  aircraft  have  been  processed 
and  reported  in  a  series  of  preceding  reports  which  are 
listed  in  Table  1.1.  However,  most  of  these  ascent  and 
descent  data,  i.e.  profile  data,  are  also  included  as  part  of 
the  general  purpose  library  described  herein.  These  profile 
data,  .summarized  in  Table  1.3,  are  contained  on  a 
separate  data  tape  from  the  radiance  arrays  and  thus  can 
be  conveniently  used  either  with  or  without  the  larger 
radiance  data  set. 

As  with  the  scanner  radiance  data,  these  profile  data 
also  have  been  organized  for  storage  by  flight  number. 
The  data  available  for  each  flight  is  summarized  in 


Table  2.3.  Due  to  a  variety  of  da'a  processing  problems, 
there  is  not  an  exact  one  to  one  correspondence  between 
the  scanner  and  profile  data  sets  listed  in  Tables  1.2  and 
1.3  The  raw  data  for  exact  correspondence  exists  on  ori¬ 
ginal  storage  tapes,  but  has  not  been  fully  retrieved  for 
inclusion  in  this  library. 

It  should  be  noted,  as  discussed  in  the  earlier  report 
series  listed  in  Table  1.1,  that  the  profile  measurements 
were  often  made  over  a  substantial  period  of  time,  and 
thus  represent  the  spatial  and  temporal  variabilities  occur¬ 
ring  along  a  flight  track  roughly  20  to  30  miles  long,  over 
a  period  of  30  to  120  minutes  in  time.  A  brief  description 
of  the  flight  profile  used  to  collect  these  data  has  been 
abstracted  from  AFGL-TR-80-0207,  Johnson  and  Gordon 
(1980)  and  is  reproduced  as  Appendix  C. 

The  radiometric  portion  of  the  profile  data  contains 
measurements  of  total  volume  scattering  coefficient  as  a 
function  of  altitude  in  each  of  the  spectral  bands  indicated 
in  Table  1.3. 

The  meteorological  portion  of  the  profile  data  con¬ 
tains  measurements  of  ambient  temperature, 
dewpoint/frostpoinl  temperature,  and  absolute  pressure. 
The  meteorological  measurements  were  made  simultane¬ 
ously  with  the  radiometric  set,  and  are  thus  listed  in  asso¬ 
ciation  with  a  specific  spectral  band  for  temporal  synchron¬ 
ization,  although  real  lime  data  at  one  second  intervals 
was  also  included  in  the  original  raw  data  sets. 

In  both  sets,  the  specific  altitude  above  ground  level 
has  been  calculated  from  measurements  of  absolute 
atmospheric  pressure,  and  the  mean  ground  elevation  as 
determined  from  local  navigation  charts 

Although  the  profile  data  points  are  all  listed  at  even 
30m  altitude  intervals,  users  should  be  aware  that  the  data 
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values  represent  the  average  of  several  measurements 
taken  throughout  that  interval.  Scattering  coefficient 
measurements  were  collected  at  a  rate  of  four  samples  per 
second,  and  meteorological  data  were  collected  at  a  rate  of 
two  samples  per  second  during  all  ascents  and  descents. 
The  aircrafts  average  rate  of  ascent/descent  throughout 
the  lower  3  km  AGL  was  approximately  1200  feet  per 
minute  (—6  meters/sec). 

Also  included  in  the  profile  data  listings  are  values 
for  several  derived  quantities  Relative  humidity  has  been 
calculated  from  the  ambient  and  dewpoint/frostpcint  tem¬ 
perature  measurements,  and  atmospheric  density  has  been 
calculated  from  the  temperature  and  pressure  measure¬ 
ments 

An  example  of  the  measured  profile  data  for  flight 
C-379  is  illustrated  in  Appendix  B.  These  plots  of  scatter¬ 
ing  coefficient,  ambient  temperature  and  computed  relative 
humidity  have  been  abstracted  from  AFGL-TR-77-0078. 
Duntley  el  a  I  (1977).  The  software  for  producing  these 
plots  is  not  included  as  part  of  this  report. 

Users  should  also  note  that  the  profile  listings  in 
this  library  do  not  include  extrapolations  above  or  below 
the  last  measured  data  point.  Thus  those  data  points 
reported  in  previous  reports  (Table  1.1)  which  were 
parenthetical  i.e.  representing  extrapolations,  are  not 
present  in  this  data  library.  Deletions  from  the  original 
data  base  for  instrumentation  reasons  are  likewise  zeroed 
out,  and  have  not  been  re-interpolated. 

THE  DATA  LISTINGS  CONTAINED  IN  THIS 
LIBRARY  HAVE  NOT  BEEN  CORRECTED  FOR 
THE  HIGH  ALTITUDE  BIAS  INDUCED  BY 
INCOMPLETE  AERODYNAMIC  PURGING  OF 
THE  AIRBORNE  INTEGRATING  NEPHELOME- 
TER  AS  DISCUSSED  IN  SECTION  2  3  OF  JOHN¬ 
SON  AND  FITCH  (1981),  AFGL-TR-82-0049 

3.0  DESCRIPTION  OF  DATA  BASE  TAPE  FORMAT 

The  primary  goals  guiding  the  development  of  this 
data  library  were  maximum  portability,  /.«?,  machine 
independence  and  simplicity  of  retrieval.  Consequently, 
both  the  tape  formal  and  the  FORTRAN  retrieval  codes 
illustrated  in  Appendix  D  and  Appendix  E  have  been 
slanted  toward  basic  simplicity  rather  than  maximum 
efficiency  in  data  packing  and  retrieval  techniques.  Both 
the  'scanner'  data  and  the  "profile"  data  are  stored  on  9 
track  tapes  at  1600  BPI  in  standard  ASCII  notation. 

As  a  cor.venience  to  the  reader,  the  title  and  intro¬ 
ductory  comment  sections  from  the  appendicized  retrieval 
routines  have  been  reoroduccd  as  Fig.  3-1  which  describes 
the  'scanner"  data,  and  Fig.  3-2  which  describes  the 
"profile'  data.  Whereas  the  information  provided  in  these 
figures  is  gene  ally  self-explanatory,  a  few  supplementary 
remarks  may  be  appropriate 


The  description  of  variable  names  in  Ftg.  3*1 
includes  several  which  may  need  clarification.  In  file  one, 
record  one,  word  S,  ’MODE",  refers  to  the  flight  pattern 
from  which  the  data  were  derived.  These  modes  are 
defined  in  Appendix  C  and  are  used  only  for  extraction 
validation  checks  at  an  earlier  level  of  processing.  The 
same  is  true  of  word  6,  "IEVENT.  This  chronological 
index  ties  the  data  array  to  a  specific  occurrence  in  the  ori¬ 
ginal  in-flight  control  log.  Neither  ’mode?  nor  'event 
numbed  information  is  necessary  for  subsequent  applica¬ 
tions  of  these  data.  Word  11,  'IFILTER",  identifies  the 
spectral  band  in  which  the  measurements  were  made  and 
is  defined  by  the  spectral  curves  shown  in  Fig.  2-1. 

In  file  two,  record  one,  word  6,  "IDENTT,  specifies 
the  beginning  of  the  measurement  interval  for  the  specific 
array  being  listed  The  times  listed  are  all  in  Greenwich 
Mean  Time  (GMT).  For  most  of  the  data  contained  in 
this  library,  the  measurement  interval  for  any  single  array 
was  180  seconds,  during  which  lime,  the  aircraft  was 
maintaining  a  straight  and  level  flight  altitude  at  an  indi¬ 
cated  airspeed  ol  150  kr.ots.  Word  9,  ’ITYPE'  defines  an 
instrumentation  set-up  which  controlled  the  airborne 
scanner's  azimuthal  rate,  and  thus  the  overall  systems 
angular  precision.  For  the  data  in  this  library,  the  type 
code  will  always  be  "3"  indicating  an  original  azimuthal 
resolution  of  0.6*  per  data  sample  which  has  been  aver¬ 
aged  to  yield  one  data  point  every  six  degrees.  Word  13, 
'1STRALIT  refers  to  whether  or  not  the  stray  light 
correction  procedure  recommended  in  an  earlier  report. 
AFGL-TR-81-027S,  Johnson  (1981)  has  been  imple¬ 
mented.  This  correction  has  not  been  made  to  the  data  in 
thi  library,  and  thus  all  arrays  should  indicate 
iSiRAL!T-0. 

In  file  two,  record  four,  there  are  two  zenith  angle 
entries  associated  with  each  data  point,  ’AVERAGE1  and 
'NOMINAL".  As  discussed  in  Johnson  (1981)  from 
which  Table  3.1  has  been  abstracted,  the  airborne 
scanner  s  operating  pattern  cal'ed  for  constant  zenith  angle 
setting,  i.e.  ’NOMINAL’  during  each  ten  second  azimu¬ 
thal  revolution  For  a  variety  of  electro-mechanical  and 
aerodynamic  reasons,  this  setting  was  not  maintained 
exactly.  The  exact  zenith  angle  analog  was  recorded  at  a  4 
sample/sec  rale  during  each  flight,  and  the  average  of 
these  actual  settings  for  each  revolution  is  listed  in  this 
data  library  as  the’  AVERAGE. 

A  similar,  bu1  shorter  series  of  comments  is  also 
appropriate  'o  Fig  3-2  which  describes  the  "profile"  data 
In  the  profile  data  comments,  file  one,  record  three,  words 
4,  5,  and  9,  MODE,  IEVENT,  and  1FLTR  have  the  same 
connotation  as  with  the  scanner  data.  However,  word  1 1 
"IH1PURGF,  refers  ic  an  upper  altitude  correction  pro¬ 
cedure  for  only  the  scattering  coefficient  data  as  discussed 
in  Johnson  (1981).  This  correction  has  not  been  made  to 
the  dau  in  this  library,  and  thus  all  profile  sets  should 
indicate  "IHIPURGF-0. 
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SCANNER  8NY/TCRRA1N  DATA  SCATTERING  COEFFICIENT  AND  RELATED 

METEOROLOGICAL  DATA 

TARE  STRUCTURE 

(1400  SRI  *  9  TRACK,  ASCII)  TAPE  STRUCTURE  t 


FILE  M  (FILE  1  GIVES  SUMMARY  OT  FLIGHTS) 

RECORO  1 t (RECORD  1  IS  RERCATEO  KFILES  TIMES) 

KFILE8 . . TOTAL  NUMSER  OF  FILES 

IYR . . . YEAR 

MONTH . MONTH 

I  DAY . DAY 

MODE . . . MODE  IDENTIFICATION 

ICVCNT . . . EVENT  NUMSER 

I  HR . . . HOUR 

MIN . .MINUTES 

ISCC . SECONDS 

ALT . .ALTITUDE  (METERS  AOL) 

IF1LTCR . FILTER  NUMSER 

FORMAT < f 19» F 10. Or  IS) 

RECORO  21 

TITLE  OF  DATA . FORMAT  (  BOA  1  ) 

RECORD  31 

FLIGHT  NUMSER i DATE* PLACE. .FORMAT (B0A1 ) 

RECORD  4! 

GEOGRAPHIC  POSITION  OF  TRACK . FORMAT < SOAt ) 

RECORD  Si 

REPORT  REFERENCES . FORMAT  <60AI  ) 

EOF 


<1400  SRIr  9  TRACK.  ASCII) 


FILE  11 

rfECORD  IS 

«tci»LSi . ro*"M<8Mn 

FLIGHT  NUMSER  r  TRACK  ID.  DATE ....  FORMAT < SOA 1 ) 
RECORD  3 1 

I  YEAR . YEAR 

MONTH . MONTH 

IDAY . . . . . DRY 

MODE . MODE  ID 

1EVENT . CVEMT  NUMSER 

I  HOUR . HOUR 

MIN . MINUTES 

1SEC . SECONDS 

1FLTR . . . FILTER  NUMBER 

NDAT . .TOTAL  NUMBER  OF  BATA  RECORBS 

IMIPUROE . HIGH  ALTITUDE  PURGE  CORRECTION 

<0«N0'1-YCB> 


FORMAT  <1119) 


FILE  21FKE  2  OIVES  UPPER  HEMISPHERE  HEADER  ONEr 
HEADER  2.A21MUTH  ANGLES r  ZCNI TH  ANGLES. 


RECORD  41 

GEOGRAPHIC  POSITION . FORMAT  <G0A1 ) 

RECORD  3! 

REFERENCE  NUMBERS . FORMAT  <BOAt) 


SCANNER 

DATA  (REPEATED  FOR  LONER 

HEMISPHERE) 

RECORD  41  (REPEAT 

RECORO  4  FOR  'NDAT'  RECORDS ) 

RECORD  I! 

JHEMI9 . . . 

DEGREES  C 

MONTH  .... 

IDAY . 

ptESS . 

I EVENT .  .  . 

(HR.HIN.SCC)OHT 
GROUND  LEVEL) 

CUBIC  METERS 

I  TYPE.  .  .  . 

(3-SCAN  RATE  OF 
AT  10  SCC/AEV) 

te  KcvoLufiOMS 

FORHAT  <12  4E 1 1  *  4  *  I 2 ) 

OREENNICH  MEAN  TIMC  (QMT) 

IAXIS . AXIS  CODE 

(2-AZIMUTH  RELATIVE  to  SUM) 

SUNAZ . SUH  AZIMUTH 

(ANOLE  IN  DECREES  TRUE) 

SUNZN . GUN  ZENITH 

(ANGLE  IN  DEGREES ) 

ISTRAL1T . SUN  ZONE  LIOHT  CORRECTION 

(0*N0. 1«YE9> 

FORMAT (A4, 41 2 r I7.F7 .0»3I2f 2f ?. 1 r 12) 


REPEAT  FILE  1  FOR  ALL  FILTERS 
EOF 


Fig.  3-2.  Scattering  coefficient  and  related  meteorological, 
format  notes. 


RECORD  2! 

FLlOHTr  (RACK  NAME 'COUNTRY r DATE . . . .FORMAT  <  01  At ) 
RECORD  3: 

1  AZIMUTHS . AZIMUTH  FROM  SUN  (40  VALUES) 

FORMAT (4014  ) 

RECORD  4: 

AVERAGE  ZENITH  ANGLE*  NOMINAL  ZENITH  ANGLE 
(IS  PAIRS  OF  VALUES) 

FORMAT ( 34F  4.2) 

RECORD  3  (REPEAT  RECORD  3*  40  TIMES) 

SCANNER  DATA . (16  VALUES) 

FORMAT ( 1SE11 *4) 

REPEAT  RECORDS  1  THRU  3  (OF  FILE  2)  FOR  SECOND  HEMISPHERE 
EOF 

REPEAT  FILE  2  FOR  KflLES'  FOR  ONE  FLIOHT 
REPEAT  FILE  1  AND  FILE  2  FOR  EACH  FLIGHT 


Tabic  3.1.  Automatic  2 if  scanner  sweep  pattern. 


Azimuth 
Revolution  No. 

COHMM 

Bevnion 

And* 

Upper  Hentopter* 
Zanlth  Angl. 

Low *r  Hemtoton 
Zaiitti  And. 

1 

2.5* 

17.5* 

92  5’ 

2 

7.5* 

12  5* 

97.5* 

3 

125* 

77.5* 

102  5* 

4 

17.5* 

72.5* 

107.5* 

5 

22.5* 

*7.5* 

112  5* 

6 

27.5* 

*2.5* 

117.5* 

7 

32  5* 

57  5* 

122.5* 

1 

37.5* 

52  5* 

127.5* 

♦ 

42.5* 

47.5* 

122.5* 

10 

47.5* 

42  5* 

137.5* 

II 

52  5* 

37  5* 

1425* 

12 

575* 

32  5' 

147.5* 

13 

27  5* 

Fig.  3-1.  Sky  and  terrain  radiance  format  notes. 
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4.0  SAMPLE  EXTRACTION  PROCEDURE 

The  dale  storage  tapes  referred  to  in  Fig.  1-1  and 
commented  upon  in  Section  3,  can  be  accessed  using  the 
computer  routines  contained  in  Appendix  D  and 
Appendix  E.  In  keeping  with  the  desire  for  simplicity,  the 
output  produced  by  these  routines  is  relatively  spartan,  but 
for  most  initial  retrievals  quite  adequate. 

An  annotated  sample  of  the  sky  radiance  array 
(UHS)  for  flight  C-378,  as  produced  by  program  READ- 
SCAN  is  illustrated  in  Fig.  4-1.  Similarly,  an  annotated 
sample  of  the  scattering  coefficient  and  related  meteoro¬ 
logical  profile  data  for  flight  C-378,  as  produced  by  pro¬ 
gram  READPROF,  is  illustrated  in  Fig.  4-2. 

4.1  Sky  A  Terrain  Radiance  Arrays 

The  following  explanatory  remarks,  directed  toward 
an  initially  naive  user,  are  related  to  the  sample  display  of 
Fig.  4-1,  and  are  complemented  by  the  format  notes  in 
Fig.  3-1. 

A.  This  introductory  table  summarizes  all  of  the  radiance 
array  data  that  is  available  for  any  given  flight.  The 
entries,  defined  in  Ale  one,  record  one,  of  Fig.  3-1, 
should  coincide  on  a  one-to-one  basis  with  the  entries 
in  Table  1.2.  Users  should  recall  that  for  every  line 
entry  in  this  introductory  table,  4- 1  A,  there  will  be  two 
data  arrays  on  the  storage  tape.  One  upper  hemisphere 
array  and  one  simultaneously  measured  lower  hemi¬ 
sphere  array. 

B.  This  two  line  array  identifies  the  sixty  discrete  azimuth 
angles  associated  with  the  radiance  data  which  follows. 
Notice  that  the  array  is  ordered  such  that  the  zero 
degrees  azimuth  is  with  reaped  to  the  solar  position, 
and  thus  for  this  array  is  equivalent  to  146  7  degrees 
true,  as  indicated  in  the  header  two  lines  above. 

C  This  two  line  array  identifies  the  eighteen  discrete  zen¬ 
ith  angles  associated  with  the  radiance  set  which  fol¬ 
lows.  Notice  that  there  are  two  entries  in  this  array  for 
each  of  the  eighteen  zenith  angle  increments  in  the 
radiance  array;  the  first  being  the  "AVERAGE1  zenith 
angle,  and  the  second  being  the  "NOMINAL".  Note 
that  the  listed  "NOMINAL"  values,  87.50,  82.50, 
77.50.  etc.  are  always  in  accordance  with  the  design 
specification  listed  in  Table  3.3,  but  for  the  reasons 
mentioned  in  Section  3,  the  true  average  values  will 
vary  slightly  during  any  given  flight. 

D.  This  two  line  array  represents  the  first  of  sixty  which 
contain  the  actual  radiance  measurements.  Each  entry 
is  identified  by  its  corresponding  azimuth  and  zenith 
angle  entry  in  B  and  C.  Thus,  since  the  radiance  data 
is  ordered  by  increasing  azimuth  angle,  this  first  group 
of  eighteen  radiance  points  are  all  at  azimuth  zero 
degrees.  The  second  group  E,  are  all  at  azimuth  six 
degrees,  etc.  as  indicated  by  the  sequence  listed  in  B. 

Within  each  azimuth  group,  it.  D  or  E,  etc.,  the 
eighteen  data  points  are  identified  in  terms  of  their 
zenith  angle  by  the  "AVERAGE'  entries  shown  in 

C. 


For  example,  consider  the  data  point  outlined  in 
group  E. 

a.  The  radiance  value  is  0.5594E+03  W/ttm1. 

b.  The  azimuth  angle  relative  to  the  sun  for  this 
data  point  is  6  degrees,  since  it  is  in  the  second 
azimuth  group  as  identified  in  B. 

c.  The  zenith  angle  for  this  data  point  is  33.95 
degrees  since  it  is  in  the  fourth  position  of  the 
second  row  as  identified  in  C.  Note  that  the 
equivalent  "NOMINAL"  value  is  27.5  degrees. 

A  default  code,  "E+23"  will  appear  as  the  last  four 
characters  in  a  radiance  value  to  indicate  that  the  radiome¬ 
ter  was  ofTscale,  beyond  its  maximum  measureabie  radi¬ 
ance. 

4.2  Scattering  Ceeftdaet  and 
Related  Meteorological  Profiles 

The  following  explanatory  remarks,  also  directed 
toward  an  initially  naive  user,  sre  related  to  the  sample 
display  of  Fig.  4-2,  and  are  complemented  by  the  format 
notes  in  Fig.  3-2. 

There  are  no  introductory  tables  associated  with  the 
profile  dau  listings  due  to  the  simplicity  of  the  listing.  For 
each  flight  there  will  be  one  profile  array  for  each  spectral 
band  in  which  scattering  coefficients  were  made,  as  listed 
in  Table  1.3.  The  appropriate  spectral  band  is  identified  in 
each  array's  header  information  as  shown  by  the  notation, 
A,  in  Fig.  4-2. 

As  illustrated  in  Fig.  4-2,  the  profile  data  is  ordered 
according  to  altitude  as  listed  in  Cel  1.  For  each  altitude 
presented,  there  will  be  an  associated  value  for  each  of  the 
parameters  shown  in  Section  F  of  Fig.  3-2,  including  the 
time  at  which  the  measurements  were  made. 

Thus,  for  example,  consider  the  row  of  data  points 
identified  by  the  notation,  B. 

From  Col.  I:  the  attitude  of  the  measurements  was  1 500m 
AGL. 

From  Col.  If:  the  measured  value  of  ambient  temperature 
was  0.7"C. 

From  Col.  IV:  the  computed  value  of  relative  humidity 
was  92%. 

From  Col.  VII:  the  measured  value  of  scattering 
coefficient  in  filter  2  was  26695E-04  m~x,  which  is  reliable 
to  only  three  significant  figures.  it.  2.67E-04  m-*. 

From  Col.  VIII:  the  time  at  which  the  measurements  were 
made  was  9  hrs.  59  min.  38  sec.  Greenwich  Mean  Time. 

Users  should  be  aware  that  although  the  values  con¬ 
tained  in  these  listings  are  recorded  to  five  significant 
figures,  as  a  formatting  convenience,  this  precision  is  far 
beyond  that  appropriate  for  the  transducers  involved. 
Experimental  reliability  beyond  two  or  three  significant 
figures  is  not  to  be  expected 
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SCATIlflNG  C'JlK'CltNT  A«('  ftLAttt'  rt(  1  COROlG&I  CAL  I'AU  -  ■  —  tel  MCI 

Flight  c-J’e  taken  at  aoditmum  t*ao  •  i»t  mbac  aom  i:  bat  1*7*  —  tel  tel 

7 A  5  12  7  11  0  54  It  MR  32  0  (A)  —tel  te) 
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1170  3.5721(400  2.8S71E-01  7.0048(401  8.7370(402  1.1024(400  2.0140(-0«  03344 
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fig.  4-2.  Annouied  lam  pie  output:  Scattering  eoeOaenl  A  retetea  meteoroioficel  data. 


5.  ACKNOWLEDGEMENTS 

Thi*  report  hu  been  prepared  for  the  Air  Force 
Geophysics  Laboratory  under  Contract  No.  F1962I-82-C- 
0060.  The  authors  wish  to  thank  the  members  of  the 
Visibility  Laboratory  technical  staff  for  their  asaistance  in 
preparing  these  data,  and  in  particular  to  acknowledge  the 
contributions  of  Ms.  Alicia  G.  Hill,  Mr.  John  C.  Brown, 
and  Mr.  James  Rodriguez,  our  specialists  in  computer 
assisted  document  preparation. 

6.0  REFERENCES 

Duntley.  S.Q.,  R.W.  Johnson.  J.l.  Oordon,  and  A.R. 
Boileau  (1970),  “Airborne  Measurements  of  Optical 
Atmospheric  Properties  at  Night,”  SIO  Ref.  70*7, 
AFCRL-70-0137,  NTIS  No.  AD  170  734,  University 
of  California,  San  Diego,  Scrippa  Institution  cf 
Oceanography,  Visibility  Laboratory. 

Duntley.  S.Q.,  R.W.  Johnson,  and  J.l.  Gordon  (1972s), 
"Airborne  Measurements  of  Optical  Atmospheric 
Properties  in  Southern  Germany,"  SIO  Ref.  72-64, 
AFCRL-72-0255,  NTIS  No.  ADA  747  490,  Univer¬ 
sity  of  California,  San  Diego,  Scripps  Institution  of 
Oceanography,  Visibility  Laboratory. 

Duntley,  S.Q.,  R.W.  Johnson,  and  J.l.  Oordon  (1972b), 
"Airborne  and  Ground-Based  Measurements  of  Opt¬ 
ical  Atmoapheric  Properties  in  Central  New  Mex¬ 
ico,”  SIO  Ref.  72-71,  AFCRL-72-0461,  NTIS  No, 


ADA  7S1  936,  University  of  California,  Sen  Diego, 
Scrippa  Institution  of  Oceanography,  Visibility 
Laboratory. 

Duntley,  S.Q.,  R.W.  Johnson,  and  J.I.  Gordon  (1975), 
“Airborne  Measurements  of  Optical  Atmoapheric 
Properties,  Summary  and  Review  II,"  SIO  Ref.  75- 
26  ,  AFCRL-TR-75-0457,  NTIS  No.  ADA  022  675, 
University  of  California,  San  Diego,  Scrippa  Institu¬ 
tion  of  Oceanography,  Visibility  Laboratory. 


Duntley,  S.Q.,  R.W.  Johnson,  and  J.l.  Gordon  (1977), 
“Airborne  Measurements  of  Atmoapheric  Volume 
Scattering  Coefficients  in  Northern  Europe,  Spring 
1976,"  SIO  Ref.  77-S,  AFGL-TR-77-0078,  OTIS  No. 
ADA  046  290,  University  of  California,  San  Diego, 
Scrippa  Institution  of  Oceanography,  Visibility 
Laboratory. 

Haring,  W.S.  0982),  "An  Operational  Technique  For 
Estimating  Visible  Spectrum  Contrast  Transmit¬ 
tance,"  SIO  Ref.  12-1,  AFOL-TR-il-0191,  NTIS 
No.  ADA  111  823,  University  of  California,  San 
Diego,  Scripps  Institution  of  Oceanography,  Visibility 
Laboratory. 

Johnson,  F.S  (1954),  "The  Solar  Constant,”  J.  Motor., 
II.  pp.  431-439. 


-10- 


likTs. 


Johnson,  R.W.  and  J.l.  Gordon  (1980),  'Airborne  Meas¬ 
urements  of  Atmospheric  Volume  Scattering 
Coefficients  in  Northern  Europe,  Summer,  1978,” 
SIO  Ref  80-20,  AFGL -TR-80-0207,  NT1S  No  ADA 
097  134,  University  of  California,  San  Diego,  Scripps 
Institution  of  Oceanography,  Visibility  Laboratory. 

Johnson,  R.W.  and  BW  Fitch  (1981),  “A  Review  of 
Measured  Atmospheric  Optical  Properties  and  Their 
Contemporary  Aerosol  Size  Distributions,”  SIO  Ref. 
82-22.  AFGL-TR-82-0049,  NT1S  No.  ADA  123  940. 
University  of  California,  San  Diego,  Scripps  Institu¬ 
tion  of  Oceanography,  Visibility  Laboratory. 


Johnson,  R.W.  (1981),  “Airborne  Measurements  of 
European  Sky  and  Terrain  Radiances,”  SIO  Ref.  82- 
2.  AFGL-TR-8 1-0275,  NT1S  No.  ADA  114  637, 
University  of  California,  San  Diego.  Scripps  Institu¬ 
tion  of  Oceanography,  Visibility  Laboratory. 


Johnson.  R.W.  and  WS  Hering  (1982).  “An  Analysis  of 
Natural  Variations  In  European  Sky  and  Terrain 
Radiance  Measurements,"  SIO  Ref.  82-6,  AFGL- 
TR-8 1-03 17,  NT1S  No.  ADA  120  487,  University  of 
California,  San  Diego,  Scripps  Institution  of 
Oceanography.  Visibility  Laboratory. 


-11- 


RADIANCE  IWft-'  *m'')  RADIANCE  (Wft -*  m‘J  (im-1) 


APPENDIX  A 


n«.  a  i. 


Flight  C-379  (Rod by,  Denmark),  dear  iky  over  ocean  with  glitter,  300m  AGL. 
(Fig  3-7  from  APOL-TK  S1-03I7) 
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F)f-  B-l.  Scattering  coefficient  A  related  meteorological  data  for  (tight  C-379. 
(Abu reeled  from  lectiona  6*7  in  AFGL-TR-77-007S) 
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FLIGHT  PROFILE  SUMMARY 
AIRBORNE  SYSTEM 

The  data  collection  sequence  for  the  airborne  system 
was  broken  into  five  standardized  elements:  (1)  preflight 
warmup  and  calibration  check,  (2)  straight  and  level 
sequences,  (3)  vertical  profile  sequences,  (4)  in-flight  cali¬ 
bration  checks,  and  (S)  post-flight  calibration  check. 

An  illustration  of  the  typical  flight  pattern,  which 
was  used  for  most  flights,  is  shown  in  Fig.  C-l.  In  this 
stylized  pattern,  two  basic  elements,  the  straight  and  level 
(ST&LV)  and  the  vertical  profile  (V-PRO),  were  com¬ 
bined  to  yield  the  tout  mission  flight  plan.  A  description 
of  these  two  pattern  elements  and  the  calibration  elements 
is  detailed  in  AFCRL-72-0255,  Duntley  era/.  (1972a), 
modiAed  in  AFCRL-TR-75-0457,  Duntley  etal.  (1975b), 
and  summarized  in  the  following  paragraphs. 

1.  Straight  and  Level  runs  (ST&LV),  Mode  03  -  The 
ST&LV  runs  were  primarily  2*r  scanner  runs.  The 
measurement  of  upper  and  lower  hemisphere  radiance 
distributions  had  top  priority.  One  sky  mode  scanner 
pattern  (192  seconds)  plus  one  sun  mode  scanner  pat¬ 
tern  (64  seconds)  were  run  at  each  altitude  with  each 
of  two  optical  filters. 

’Sky  Mod^  and  ’Sun  Mode/1  refer  to  procedural  adjust¬ 
ments  which  were  made  to  enable  radiometric  meas¬ 
urements  of  acceptable  fidelity  in  the  high  gradient 
regions  of  the  solar  aureole.  In  *Sky  Mode' 
configuration,  radiometric  sensitivity  was  adjusted  to 
optimize  measurement  of  sky  radiances  far  removed, 
e.g.  >30*.  from  the  solar  disc.  Conversely,  in  ’Sun 
Modi',  sensitivity  was  reduced  to  enable  direct  solar 
disc  and  aureole  measurements  at  the  expense  of  the 
off-sun  regions. 

During  ST&LV  runs  the  aircraft  maintained  a  fixed 
heading,  a  constant  indicated  airspeed  of  1 50  knots  or 
less,  and  a  2.5  degree  nose-high  flight  altitude. 

2.  Vertical  Profile  runs  (V-PRO),  Mode  07  -  The  V-PRO 
runs  were  primarily  integrating  nephelometer  and  vari¬ 
able  path  function  meter  runs.  The  measurement  of 
the  total  scattering  coefficient  profile  had  top  priority. 
Second  priority  was  measurement  of  the  vertical  path 
function  profile.  Each  V-PRO  ascent  or  descent  was 
made  using  a  single  filter. 

During  the  V-PRO  runs  the  aircraft  maintained  a  Axed 
heading,  with  the  sun  off  the  left  wingtip,  and  a  flight 
attitude  not  exceeding  4  degrees  nose  down  or  8 
degrees  nose  up.  An  average  rate  of  climb  or  descent 
of  1200  feet/minute  was  optimum,  and  airspeed  was 
not  critical,  but  remained  constant  once  established. 

3.  Cross-Calibration  Climbs  (X-CAL),  Mode  08  -  The 
X-CAL  climbs  were  specifically  designed  to  validate  the 
performance  of  the  UHS,  LHS,  and  ERT  radiometer 


systems.  The  simultaneous  measurement  of  a  com¬ 
mon  uniform  segment  of  sky  by  these  three  radiome¬ 
ters  had  top  priority.  Two  X-CAL  climbs  were  associ¬ 
ated  with  each  standard  profile,  one  preceding  the  first 
ST&LV  run  and  the  second  following  the  last  ST&LV 
run.  Both  sky  mode  and  sun  mode  measurements 
were  made  with  the  UHS  system. 

During  the  4-minute  X-CAL  climb  the  aircraft  main¬ 
tained  a  fixed  heading,  with  the  sun  in  the  aft  hemi¬ 
sphere,  and  a  5-degree  nose-high  flight  attitude.  The 
aircraft  was  flown  directly  toward  the  dearest  and  most 
uniform  portion  of  the  sib'  as  waa  practical. 

4.  Calibration  Blocks  (A/D  CAL),  Mode  00,  M-CAL, 
Mode  01,  N-CAL,  Mode  09  -  The  32-aecond  blocks  of 
calibration  data  were  inserted  periodically  throughout 
the  entire  data  mission.  They  were  designed  to  provide 
calibration  update  information  to  the  poet-flight  com¬ 
puter  processing  sequences.  There  are  21  assorted  cali¬ 
bration  Mocks  associated  with  each  standard  profile. 

During  these  calibration  blocks  there  were  no  project- 
imposed  requirements  upon  aircraft  speed  or  attitude. 

GENERAL  FLIGHT  PATTERN 

In  this  profile,  ST&LV  daU  runs  were  made  using 
two  different  spectral  filters  at  each  of  four  altitudes.  The 
ascent  V-PRO  waa  made  using  the  first  of  the  two  Alters, 
and  the  descent  V-PRO  waa  made  using  the  second.  After 
the  descent  V-PRO,  the  entire  sequence  was  repeated 
using  a  second  pair  of  filters. 

The  idealized  flight  profile  would  result  in  all  ground 
tracks  falling  on  a  single  line  running  between  the  Initial 
Point  (I.P.)  and  the  Turning  Point  (T.P.).  In  practice,  the 
ST&LV  elements  were  actually  stacked  in  a  slab  of  atmo¬ 
sphere  approximately  30  miles  (48  km)  long,  0.5  mile  (0.8 
km)  wide,  and  4  miles  (64  km)  high. 

Periodically,  in  response  to  specialized  data  require¬ 
ments  or  weather  conditions,  supplementary  flight  patterns 
were  added  to  the  mission  profile.  For  example,  a  pattern 
made  up  of  a  (2+3)  profile,  /.e.,  two  spectral  filters  at 
each  of  three  altitudes  waa  used  as  was  a  (2+2)  profile, 
l.e. ,  two  spectral  filters  at  each  of  two  altitudes.  Both  the 
(2+3)  and  (2+2)  profiles  are  generally  considered  low  to 
medium  altitude  profiles,  and  were  normally  used  on 
flights  performed  under  a  full  overcast  or  low  to  inter¬ 
mediate  level  cloud  decks.  Three  flights  consisted  of 
V-PRO  climb  outs,  supplemented  with  only  directional 
scattering  measurements  at  the  minimum  and  maximum 
altitudes. 

At  the  conclusion  of  each  mission,  the  radiometric 
data  which  were  recorded  and  stored  on  magnetic  tape 
were  returned  to  the  Visibility  Laboratory  for  computer 
reduction  and  analysis 
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C* ******************************************************** 

c  title: 

C  READSCAN... . FORTRAN  PROGRAM 

r:  ***************************** **************************** 

c  purpose: 

C  READ  SCANNER  TRANSMITTAL  TAPE  OF  SUN/TERRAIN  DATA 

C****************** *************************************** 


SCANNER  SKY/TERRAIN  DATA 


TAPE  STRUCTURE 
<1600  DP  I r  V  TRACK  >  ASCII) 
miLL  RECORDS  240  DYTES)  • 


NUMBER 


TIMES) 
OF  FILES 


FILL  It  -FILE  1  GIVES  SUMMARY  OF  FLIGHTS) 

RECORD  It (RECORD  1  IS  REPEATED  ' KF ILES '  TIMES) 

KTILEC .  . TOTAL  NUMBER  OF  FILES 

1YR . YEAR 

MONTH  .  . .  .MONTH 

I  DAY . DAY 

MODE . MODE  IDENTIFICATION 

I  EVENT . EVENT  NUMBER 

I  HR . HOUR 

MIN  .  . MINUTES 

I  SEC . SECONDS 

ALT . ALTITUDE  (METERS  AGL  ) 

IFILTER . FILTER  NUMBER 

FORMA T(?I5fF10.0rI5) 


RECORD  2*. 

TITLE  OF  DATA . FORMAT  (  B0A1  ) 

RECORD  3J 

FLIGHT  NUMBER  »  DATE »  PLACE . . FORMAT ( 80A1 ) 

RECORD  4J 

GEOGRAPHIC  POSITION  OF  TRACK . FORMAT ( 80A1 ) 

RECORD  5: 

REPORT  REFERENCES . FORMAT  (80 A 1  ) 


FILE 


FILE  2  GIVES  UPPER  HEMISPHERE  HEADER  ONE. 
HEADER  2  >  AZ I MUTH  ANGLES  t  ZENITH  ANGLES  r 
SCANNER  DATA  (REPEATED  FOR  LOWER  HEMISPHERE) 


RECORD  i: 

JHEMIS . .  .  .HEMISPHERE 

IYR . YEAR 

MONTH . MONTH 

I  DAY . .  .DAY 

I  EVENT . EVENT  NUMBER 

IDENT2 . DATA  START  TIME  (HR, MIN. SEC 

ALT . ALTITUDE  (METERS  ABOVE  GROUND  LEV 

IFLTR . FILTER  NUMBER 

I  TYPE . TYPE  CODE 

( 3=SCAN  RATE  OF  18  REVOLUTI 
AT  10  SEC/REV) 

I  AX  I . . AXIS  CODE 

( 2*AZ I MUTH  RELATIVE  TO  SUN) 


(  HR , MI N , SEC )GHT 
GROUND  LEVEL) 


REVOLUTIONS 


V.  A 

»  •  *  .  'V 

•*.nV»V 


?::W 


v- \>y 


•  * .  •  *  * 

PC 


-;.v. 


•  * . 
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C 

c 

c 

c 

c 

c 

c 

L 


SUNAZ . SUN  AZIMUTH 

< ANGLE  IN  DEGREES  TRUE) 

SUNZN . SUN  ZENITH 

(ANGLE  IN  DEGREES) 

ISTRALIT . SUN  ZONE  LIGHT  CORRECTION 

( 0=N0 1 1 =YES  > 


FORMAT (A4» 412. 17, F7. 0.312. 2F7. 1 » 12) 
RECORD  2: 


C  FLIGHT  » 1  RACK  NAME , COUNTRY , HATE  ....  FORMA T ( 8 1 A 1 ) 

C  RECORD  3t 

C  IAZtMUTHS . AZIMUTH  FROM  SUN  (60  VALUES) 


C 


FORMAT ( 6014 ) 


C  RECORD  Al 

C  AVERAGE  ZENITH  ANGLE.  NOMINAL  ZENITH  ANGLE 

C  (18  PAIRS  OF  VALUES) 


C 

C 

c 

c 

»: 


FORMAT ( 36F6 . 2  > 

RECORD  5  (REPEAT  RECORD  5.  60  TIMES) 

SCANNER  DATA . (18  VALUES) 

FORMAT  < 13E 1 1 .4) 


C  REPEAT  RECORDS  1  THRU  S  < OF  FILE  2)  FOR  SECOND  HEMISPHERE 

C 

U  EOF 

C 

C  REPEAT  FILE  2  FOR  ' KFILES '  FOR  ONE  FLIGHT 

C 

C  REPEAT  FILE  1  AND  FILE  2  FOR  EACH  FLIGHT  T* 

C* ******************************************************** 

C  PROGRAMMER 

C  MIRIAM  K,  OLEINIK 

C  VISIBILITY  LABORATORY 

C  SCRIFPS  INSTITUTION  OF  OCEANOGRAPHY 

C  UNIVERSITY  OF  CALIFORNIA  SAN  DIEGO 

C* ************************************* ********** ********* 

C  ****  NOTE  **** 

C  TAPE  INPUT  FUNCTIONS  ARE  COMPUTER  SYSTEMS  DEPENDENT 

C 

L****************** ****************  ************************ 
DIMENSION  ZA( 18) »ZN( 18) . I AZ IMUTH ( 60 ) . I T I TLE ( 240 ) . 

1  DATA  (18).  JERROR <  240 ) 

DATA  NFLIGHT/1/ 

DATA  LUU/6/ 

C 

C  NFLIGHT  =  TOTAL  FLIGHTS  TO  READ 

C  LUW  "  LOGICAL  UNIT  FOR  OUTPUT  PRINTER 

C 


C  LUTAPl  IS  LOGICAL  UNIT  OF  ID  RECORDS 

C  LUTAPE  IS  LOOICAL  UNIT  OF  DATA  RECORDS 

C 

LUTAP1=8Q 
LUT  APE  =  60 

C  LOOP  FOR  TOTAL  FLIGHTS 

C 

=  1. NFLIGHT 


DO  4000  IFLIGHT 
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C 

C  OPEN  FILE  FOR  ID  RECORDS 

r 

OPEN<UNIT=LUTAPl * RECORDSI ZE=240 * TYPE= ' OLD '  ) 

C 

C  READ  IDENTIFICATION  RECORD  <1> 

C 

IFILE  =  1 

100  FORMAT  <9I5*F10,0*I5) 

200  CONTINUE 

READ<LUTAP1 1 100  * ERR  =  320 ) KF ILES *  I YR *  MONTH  * I  DAY .MODE. 

1  IEVENT, IHR*MIN, ISEC* ALT* I FILTER 
IFILE  =  IFILE  +  1 

URI TE (LUU* 300 JKFILES* I YR* MONTH* I  DAY » MODE * IEVENT *IHR* 
1  MIN  * ISEC  *  ALT  *  IF ILTER 
300  FORMAT  (1X*9I5*F10.0*I5) 

GO  TO  340 

C 

C  ERROR  IN  READ *  READ  IN  ALPHA  AND  PRINT 

C 

320  URI TE ( LUW *  7010 ) 

BACKSPACE  LUTAP1 

READ(LUTAP1 *970) ( JERROR( JA) *-A=l *60) 

UR  I TE ( LUU  t 550 ) <  JERROR ( JA  > »JA  =  1 »60) 

C 

C  TEST  IF  END  OF  ID  RECORD  (1)  READ 

C 

340  IF< IFILE. GT.KFILESIQO  TO  400 
C  GO  READ  NEXT  ID  RECORD 

GO  TO  200 
400  CONTINUE 

C 

C  READ  IDENTIFICATION  RECORD  (2) 

C 

DO  600  IHEAD  =  1>4 

READ ( LUTAP 1 » 500) ( IT ITLE < JA ) t JA= 1 r 80 ) 

500  FORMAT  <  240A1 ) 

550  FORMAT (3< lX»80Al/> ) 

600  UR  I  TE  <  LUU » 550  )  ( I T I TLE  (  JA ) »  JA=*  1 »  80 ) 

C 

C  LOOP  FOR  TOTAL  FILES  THIS  FLIGHT 

C 

DO  3000  IFILE  =  1»  KFILES 
C 

C  OPEN  FILE  FOR  DATA  RECORDS 

C 

OPEN  ( UNI  T  =  LUTAPE  » RECORDS  I  ZE«240»  TYPE«=' OLD'  ) 

C 

C  LOOP  FOR  2  HEMISPHERES  (UPPER  AND  LOWER) 

C 

DO  2000  IHEMIS  =  1*2 


C 

C 

C 


READ  SCANNER  HEADER  <1> 


CJOCJ  (j  (j  U  U  J  <J  u  c  U  CJCJCJ 
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READ ( LUTAPE *  700 tERft-740) JHEKIS. I YR» MONTH » IDAY. IEVENT » 
1  IDENf2.ALT»IFLTR.ITYPE»IAXIS»  SUNAZ  » SUNZN  . ISTRALIT 

700  FORMAT (1X.A3. 4I2.I7.F7. 0.312. 2F 7.1 » 12) 

WRITE(LUW»700)JHEMIS.IYR»  MONTH  .IDAY.IEVENT.IDENT2. 

1  ALT. IFLTR. I  TYPE.  I  AXIS. SUNAZ » SUNZN » ISTRALIT 
GO  TO  750 
C 

C  ERROR  IN  TAPE  READ.  READ  IN  ALPHA  AND  PRINT 

C 

740  WRITE(LUW.7010) 

BACKSPACE  LUTAPE 

READ (LUTAPE .970) ' JERROR ( JA )» JA= 1 . 51 > 

WRITE(LUW»  550  >( JERROR ( JA ) » JA=1»51 ) 

READ  SCANNER  HEADER  <2> 

750  REALi<  LUTAPE  .500)  <  ITITLE<  JA)  .  JA=1  »  81  ) 

WRITE (LUU. 500)  ( I  TITLE ( JA) . JA=1 .81  ) 

AZIMUTHS 

RE AD (LUTAPE. 800. ERR=845 ) < IAZIMUTH( JA) . JA=1 .60) 

800  F  ORMAT ( 60  14) 

URITE(LUU,840) ( I AZ I MUTH < JA ) . JA= 1 . 60 ) 

8 4 C  FORMAT (2(30I4/>  > 

GO  TO  850 

ERROR  IN  TAPE  INPUT.  READ  IN  ALPHA  AND  PRINT 

845  WRITE (LUU. 7010) 

BACKSPACE  LUTAPE 

RE AD (LUTAPE. 970) ( JERROR ( JA >. JA= 1 . 240 > 

WR1TE(LUW.550) < JERROR ( JA ). JA" 1 . 240 ) 

ZENITHS 

850  READ (LUTAPE. 855. ERR  =  B65 ) <ZA< JA) »ZN<  JA) , JA*1. 18) 

855  FORMAT ( 36F6 • 2 ' 

WRITE (LUU. 860) ( Z A ( JA > . ZN < J A ) , JA=1 . 18) 

860  FORMAT  <2(18F7.2/>> 

GO  TO  890 

865  WRITE (LUW. 7010) 

BACKSPACE  LUTAPE 

READ (LUTAPE .970)  (JERROR ( JA) . JA=1 .216) 

WRITE (LUW. 550) ( JERROR ( JA > . JA= 1 . 2 1 6 ) 

890  CONTINUE 

LOOP  FOR  SCANNER  DATA 

DO  1000  IDATA  =  1.60 

READ ( LUTAPE. 900. ERR  =  960 X  DATA ( JA). JA= 1.18) 

900  FORMAT ( 1 8E 1 1 . 4  ) 

WRITE(LUW.910) (DATA( JA) . JA=1 , 18) 

910  FORMAT (2( 1X.9E11 . 4/) ) 

GO  TO  1000 
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ERROR  IN  TAPE  READ »  READ  IN  ALPHA  AND  PRINT 

960  WRITE(LUUt7010) 

BACKSPACE  LUTAPE 

READ (LUT APE » 970) <  JERROR < JA >» JA= 1 » 1 9B ) 

970  FORMAT ( 240A1 ) 

WRITE (LUW»S50) (JERROR ( JA) » JA= It  198) 

1000  CONTINUE 

GO  READ  LOWER  HEMISPHERE 

2000  CONTINUE 

CLOSE  AND  SAVE  DATA  FILE 

CLOSE < UNI T=LUT APE. DISPOSE*' SAVE' > 

LUTAPE=LUTAPEF1 

GO  READ  NEXT  FILE  THIS  FLIGHT 

3000  CONTINUE 

CLOSE  AND  SAVE  HEADER  RECORDS 

CLOSE ( UN IT  =  LUT API » DISPOSE* ' SAVE  '  ) 

LUTAP1-LUTAP1F1 

GO  READ  NEXT  FLIGHT 

4000  CONTINUE 

ERROR  BRANCH  MESSAGE 

7010  FORMAT <  IX »  '  ****  ERROR  IN  TAPE  READ f  READ  IN  ALPHA  AND  PRINT  ') 
9000  CONTINUE 
END 
I 
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CtSttHtf «**¥**« *********** *********** *«*>****¥«¥*** 

TITLE 

READPROF .FOR 


PURPOSE 

READ  SCATTERING  COEFFICIENT  AND  RELATED 
METEOROLOGICAL  DATA 
(TRANSMITTAL  TAPE  DATA) 


C***************************************************** 

C  SCATTERING  COEFFICIENT  AND  RELATED 

C  METEOROLOGICAL  DATA 

C 

C  TAPE  STRUCTURE { 

C 


C  (1600  BPI.  9  TRACK »  ASCII) 

C  (ALL  REC0RDS=80  BYTES) 

C 

C  FILE  i: 

C  RECORD  1 t 

C  TITLE .  FORMAT  (  80A1 ) 

C  RECORD  2: 

C  FLIGHT  NUMBER *  TRACK  ID*  DATE . . . . FORMAT ( 80 A 1 ) 

C  RECORD  3: 


C  I  YEAR . YEAR 

C  MONTH . MONTH 

C  IDAY . DAY 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


MODE .  i  p  .  * . .  i  >  •  >  i  1 1  •  i  i  i  •  i  •  i  •  i  >  i  i  MODE  I D 

IEVENT »»•••••*♦•**••»»••«•••••*»•*• . EVENT  NUMBER 

IHOUR. . . <HOUR 

MIN......... . . . .  ...  MINUTES 

ISEC.. ........a.... »«.••»•.••.... ...SECONDS 

IFLTR . ••••..<••••••<••••••<••<•• .FILTER  NUMBER 

NDAT . . . TOTAL  NUMBER  OF  DATA  RECORDS 

IHIPURGE . HIGH  ALTITUDE  PURGE  CORRECTION 


( 0  =  N0 1 1  =  YES ) 


FORMAT( 1 115) 


C  RECORD  4J 

C  GEOGRAPHIC  POSITION . FORMAT  (80A1) 

C  RECORD  55 


C  REFERENCE  NUMBERS . FORMAT  (80A1) 

C  RECORD  65  (REPEAT  RECORD  6  FOR  'NDAT'  RECORDS) 

C  I  ALT . ALTITUDE  IN  METERS  (AGL) 

C  TEMP . TEMPERATURE  IN  DEGREES  C 


C  DEWTEHP 

C 

C  R  H  .  *  .  .  . 

C  PRESS.. 

C  DENS... 

C 

C  SCAT**. 

C 


DEWPOINT/FROSTPOINT  IN 
DEGREES  C 

RELATIVE  HUMIDITY (PERCENT ) 
PRESSURE  IN  MILLIBARS 
DENSITY  IN  KILOGRAMS  PER 
CUBIC  METERS 

SCATTERING  COEFFICIENT  IN 
INVERSE  METERS 


C 


I T I  ME 


TIME  (HR. MIN. SEC.  ) 


C 


GREENWICH  MEAN  TIME  (GMT) 
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C  F0RMAT(I5»6E11.4»I7) 

C 

C  REPEAT  FILE  1  FOR  ALL  FILTERS 

C  EOF 

C  PROGRAMMER 

C  MIRIAM  K.  OLEINIK 

C  VISIBILITY  LABORATORY 

C  SCRIPPS  INSTITUTION  OF  OCEANOGRAPHY 

C  UNIVERSITY  OF  CALIFORNIA  SAN  DIEOO 

C  ****  NOTE  **tr 

C  TAPE  INPUT  FUNCTIONS  ARE  COMPUTER  SYSTEMS  DEPENDENT 

c««****************t********«t*«»t**t**»tt**t**«*stttt**t* 

c 

DIMENSION  DATA(6).ITITL£<80) 

DIMENSION  IPOS (80) f IREF ( 80 > 

DATA  LUW/6/ 

C 

C  LUW  -  LOGICAL  UNIT  FOR  PRINTER 

C  LUTAPE  *  LOGICAL  UNIT  FOR  INPUT  FROM  TAPE 

LUTAPE  o  60 
C 

40  FORMAT (80AX ) 

C  OPEN  INPUT  DISK  FILE 

C 

OPENCUNIT = LUTAPE t RECORDS I ZES 80 r  TYPE* ' OLD ' ) 

C 

C 

C  LOOP  FOR  ALL  FILTERS 

C 

60  CONTINUE 
C 

C  READ  PROFILE  DATA 

C 

C  INPUT  TAPE  HEADER  INFORMATION 

C 

READ (LUTAPE* 40 » END-3000 ) < IT  I TLE  <  JA ) ► JA-1 .80) 

WRITE ( LUW » 40 > (ITITLE< JA> .JA-1 .80) 

READ (LUTAPE. 40. END-3000) < I  TITLE! JA) » JA-1. 80) 

WRITE (LUW .40) ( IT  I TLE < JA) .  JA-1. 80) 

READ ( LUTAPE » 100. END-3000) I YEAR. MONTH. I  DAY .MODE. IEVENT. 
1  I HOUR .MIN . I SEC .IFLTR.NDAT. IHIPURGE 
100  FORMAT ( 1 1 IS ) 

WR I TE( LUW. 100) I YEAR. MONTH. IDAY. MODE. IE VENT. 

1  IHOUR. MIN. ISEC. IFLTR.NDAT. IHIPURGE 
READ ( LUT APE. 40. END- 3000 ) (IPOS ( JA). JA-1 .80) 

READ (LUTAPE .40 .END-3000) ( IREF( JA> . JA-1 .80) 

WRITE (LUW .40) (IPOS ( JA) » JA-1. 80) 

WRITE (LUW. 40) (IREF ( JA) » JA-1 .80) 

DO  2000  IDATA  -  l.NDAT 

READ ( LUTAPE. 400 . END- 3000 ) I ALT . ( DATA ( JA ) »  JA-1 .6) . I  TIME 
WRITE (LUW. 400) I ALT. (DATA< JA) . JA-1 . 6 )» I  TIME 
400  FORMAT (1X.I4.1P6E11«4.I7) 

2000  CONTINUE 
C 
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READ (LUTAPE*  400 »  END® 3000) I ALT  t  < DATA  ( JA) » JA=1 t  6 ) » I  TIME 
WRITE  <  LUW 1 400 ) I ALT  » ( DATA ( JA) »  JA*1 » & ) r IT IME 
400  FORMAT <1X>I4>1P£E11 • 4  > I  7 } 

2000  CONTINUE 

00  FOR  NEXT  FILTER 

00  TO  60 
3000  CONTINUE 

WRITE (LUW  » 4000  > 

4000  FORMAT  <  IX  r '  END  OF  FILE  ') 

CLOSE  INPUT  FILE 

CLOSE<UNIT  ■  LUTAPEf DISPOSE® ' SAVE ' ) 

END 
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